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Abstract. We perform randomized benchmarking on neutral atomic quantum bits 
(qubits) confined in an optical lattice. Single qubit gates are implemented using 
microwaves, resulting in a measured error per randomized computational gate of 
1.4(1) x 10 -4 that is dominated by the system T2 relaxation time. The results 
demonstrate the robustness of the system, and its viability for more advanced quantum 
information protocols. 
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1. Introduction 

Quantum information processing has the potential to revolutionize computation and 
communication. However, the physical realization of a quantum information processor 
places stringent demands on the efficiency and reliability of the individual quantum 
gate operations that drive the computation pQ. Although many methods have been 
developed to characterize the performance of quantum gates, these may be difficult 
or computationally intensive to implement, and may be unable to distinguish between 
initialization, read-out, and gate errors 0,[3]. 

Randomized benchmarking has been proposed as a way to efficiently evaluate the 
performance of a quantum gate [4j. The basic protocol is to apply a random sequence 
of unitary operations, and measure the decay in the average output state fidelity as a 
function of the number of applied gates. This procedure has a number of advantages over 
other evaluation methods. The randomization of the sequence depolarizes the noise in 
the system, and ensures that the aggregate result is independent of a particular gate or 
proper subset of gates. Moreover, randomized benchmarking distinguishes between the 
error of the quantum gate and errors due to initialization and state detection. Finally, 
this protocol determines the average fidelity of the gate at an arbitrary point in a 
sequence of operations, establishing independence of the error rate to changes in the 
experiment during the course of the quantum computation. Randomized benchmarking 
has recently been used to evaluate quantum gates in trapped atomic ion jU [5] , liquid- 
state NMR [6j, and solid-state qubits 0IE]. 

We implement randomized benchmarking on ultracold atomic qubits trapped 
in an optical lattice. The single qubit gates are driven by microwave pulses with 
precisely controlled frequency, duration and phase. The measured error per randomized 
computational gate of 1.4(1) x 10 -4 is dominated by decoherence of the trapped atomic 
spins, and is indicative of excellent control over the applied gate operations. We expect 
that straight-forward upgrades to the system will further improve both the coherence 
of the qubits and the fidelity of the quantum gates. 

2. Experimental Setup 

We prepare an ultracold cloud of rubidium atoms ( 87 Rb) by first passing atoms from 
a hot sample through a Zeeman slower [9], and then laser-cooling them in a six- 
beam magneto-optical trap [10]. The atoms are then confined by a Ioffe-Pritchard- 
style magnetic trap, and are cooled by forced radiofrequency (rf) evaporation [TT] . 
Subsequently, a 1.5 fim dipole- force trap and quadrupole magnetic field retain the atoms, 
where the addition of the quadrupole magnetic field assists in supporting the atoms 
against gravity and allows for confinement along the propagation direction of the dipole 
beam to be comparable to that provided by the beam in the transverse directions. Here 
the atoms are further evaporatively cooled by gradually decreasing the intensity of the 
1.5 /xm light, resulting in a spin-polarized Bose-Einstein condensate [12]. 
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Transitions between the ground level hyperfine manifolds are driven by applying 
microwave radiation at frequencies near 6.8 GHz. The required frequencies are generated 
by mixing a low frequency (0 to 180 MHz) output from either a direct-digital synthesizer 
(DDS) or an arbitrary waveform generator (AWG) with the frequency-doubled output of 
a microwave signal generatorEI The use of a DDS or AWG enables simple control of the 
frequency and phase of the radiation; here, the AWG is used for efficient frequency 
sweepscl The mixer output passes through a TTL-controlled rf switch (producing 
approximately square-shaped pulses), then is amplified and sent to a microwave horn 
directed at the trapped atoms. 

Before loading the atoms into the optical lattice, we reduce the total number of 
trapped atoms by sweeping the frequency of applied microwave radiation near 6.8 GHz 
to transfer a fraction of the atoms from the trapped 2 S r 1 / 2 |F = l,ra j F = — 1) = |1, — 1) 
state to the anti-trapped 2 Si/ 2 \F = 2,rnp = —2) = |2, —2) state. Here, F is the total 
angular momentum of an atom, and mp is its projection along a quantization axis 
defined by an external magnetic field. Using this microwave sweep technique, we are 
able to prepare an ultracold atomic sample of approximately 2 x 10 4 atoms. 

The atoms are adiabatically loaded into a three dimensional, red-detuned (A/ = 
810.5 nm) optical lattice by exponentially increasing the intensity of the lattice beams to 
a maximum value over a duration of about 200 ms. In two directions, retro-reflection of 
an intensity-stabilized beam results in a lattice with spacing between sites of A//2 = 405 
nm; the optical configuration is similar to the one described in Ref. p2]. The lattice 
in the third dimension is generated by two beams (not intensity-stabilized) propagating 
at an angle of about 160°, resulting in a lattice spacing of approximately 411 nm. The 
depth of the potential between adjacent lattice sites is measured to be about 30i? r , 
determined by pulsing the lattice beams and observing an atom diffraction pattern in 
time-of-flight [Hj. Here the recoil energy is defined by E r /h = h/(2Xfm) « 3.5 kHz, 
where h is Planck's constant and m is the mass of the confined atom. The low initial 
atom number and measured lattice depth results in no more than one atom per lattice 
site, as confirmed by microwave spectroscopy [T5j . 

3. Randomized Benchmarking Experiment 

We perform randomized benchmarking on the atoms confined in the optical lattice using 
a series of microwave pulses, following the procedure described in Ref. [4J. Our qubit 
states are the magnetic field-insensitive |2,0) and |1,0) states, defined as |0) and |1), 
respectively. Immediately after being loaded into the optical lattice, the confined atoms 
are in the state |1, — 1). A microwave 7r-pulse between |1, — 1) and |0) initializes the state 
of the qubits, as illustrated in Fig.[l](a)0 A sequence of randomized computational gates 

a Rhode & Schwartz SMT 06. 

b Agilent 33250A Function/ Arbitrary Waveform Generator. 

c The low frequency signal sent to the mixer for the state initialization and state read-out 7r-pulses is 
produced by a commercial DDS signal generator: Novatech 409A. 
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Figure 1. Randomized benchmarking procedure, (a) State initialization. Atoms 
loaded into the lattice are initially in the state |1, — 1). A microwave 7r-pulse prepares 
the atoms in the qubit state |0). (b) Randomized computational gate sequence. 
A series of random microwave 7r-pulses (Pauli randomizations, PR) and 7r/2-pulses 
(computational gates, CG) is applied between the qubit states |0) and |1). A final, 
calculated 7r/2-pulse R returns the qubit to an eigenstate of <j z for measurement, (c) 
State read-out. A microwave 7r-pulse maps |0) to |1,— 1) in order to read out the 
relative occupation of the qubit states in a Stern-Gerlach configuration. 



is then applied between the |0) and |1) qubit states (Fig. [11(b)). A single randomized 
computational gate is defined as the concatenation of a randomized 7r-pulse or identity 
operation (Pauli randomization, PR) and a 7r/2-pulse (computational gate, CG). An 
ideal Pauli randomization has the functional form e ±iap7T ^ 2 , where a p are the standard 
Pauli operators with p = l^x^y^z. The sign ± and Pauli operator a p are chosen at 
random with uniform probability at each occurrence in the sequence. As in previous 
implementations, the a z Pauli randomization pulses are implemented only by changing 
the frame for all subsequent pulses [U El [7] . Although no pulse is applied for the aj 
and a z operations, the duration of the operation in the sequence is kept the same across 
all Pauli randomizations. The form of an ideal computational gate is e ±l(7c7T ^ ) with the 
sign ± and c = x,y each chosen at random with equal probability. Although some 
previous implementations of randomized benchmarking have included the change of 
frame operation ( e ±zcr W 4 ) as a computational gate [6], here we follow Ref. pQ El El [8] 
and only implement x ) y rotations as computational gates. Following the sequence of 
randomized computational gates, a final 7r/2-pulse with random sign (nested between 
two more PR pulses) is performed, that is calculated to return the qubit state to an 
eigenstate of a z in the absence of errors. The qubit state is read out by mapping |0) to 
|1, — 1) with a microwave 7r-pulse (Fig. [1(c)), lowering the depth of the optical lattice, 
and releasing the atoms in a magnetic field gradient. This Stern-Gerlach type analysis 
spatially separates the different spin states 1) and |1) (= |1,0)) in time-of-flight, 
allowing us to determine the relative population in each by absorption imaging. 

As the number of randomized computational gates implemented in a sequence is 
increased, the accumulated gate error reduces the measured single qubit fidelity T of 
the output state. The fidelity is defined as the overlap of the ideal (pideai) and measured 
(p) density matrices, so that T = tr(pid ea ip)- The average output state fidelity (J 7 ) 
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Figure 2. Fidelity of randomized benchmarking sequences versus the number of 
randomized computational gates (CG). (a)-(d) Here, SI, S2, S3, and S4 denote 
the measurements obtained for the four different random computational sequences, 
each combined with eight Pauli randomization sequences. Averaging each over all 
Pauli randomization sequences and fitting (solid lines) to Eq. Q] results in d = 
{2.6(4), 2.9(3), 2.9(3), 2.3(3)} x 1(T 4 for {SI, S2, S3, S4}. (e) Measured output state 
fidelity for all combinations of computational and Pauli randomization sequences. The 
scatter in the data suggests that coherent errors may contribute to the decrease in 
fidelity, (f) The output state fidelity averaged over all randomized gate sequences. 
The results are fit to Eq. ffl yielding d if = 1.8(2) x 10" 2 and d = 2.7(2) x 10" 4 , for 
an error probability per randomized computational gate of £ g = d/2 = 1.4(1) x 10 -4 . 
Error bars are statistical, and represent the standard deviation of the mean. 



decays exponentially as j4] 

7 = \ + \iX-dij)ix-d) 1 (i) 

where dif is the probability of depolarization due to initialization and read-out, d is the 
average probability of depolarization of a single randomized computational gate (PR 
and CG), and I is the number of randomized computational gates applied. The average 
decrease in fidelity per randomized computational gate (average error per gate) £ g is 
related to the probability of depolarization d by [4J 
d 

£ 9 = 2> (2) 
giving a fidelity of 1/2 for a completely depolarized qubit. Thus, we can apply sequences 
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of varying lengths of randomized computational gates to our single qubits, and fit the 
decay in the average output state fidelity to Eq. [I] to determine d and S g . 

We applied a total of 32 randomized computational gate sequences, derived from 
4 random sequences of computational gates and 8 sequences of Pauli randomizations. 
Each randomized computational gate sequence was truncated at 15 different lengths {1, 
2, 3, 5, 8, 13, 21, 34, 55, 89, 145, 235, 380, 615, 995}, for a total of 4 x 8 x 15 = 480 
sequences. The measured output state fidelity for each of these 480 implementations 
is shown in Fig. E](a)-(e). The average fidelity of the 32 sequences at each truncation 
length is given in Fig. [2(f), and a weighted fit of this data to Eq. [1] yields an average 
error per gate £ g = d/2 = 1.4(1) x 10" 4 . 

4. Evaluating Errors 

We performed several experiments to determine whether control errors or decoherence 
limit the measured error per computational gate. The two primary control parameters 
are the frequency and duration of the applied microwave pulses. The frequency to drive 
transitions between the two qubit states is generated by mixing the (high frequency) 
signal from a microwave signal generator with the (low frequency) output of a custom 
DDS boardH The frequency and phase of the DDS are digitally controlled by a custom 
field-programmable gate array (FPGA) boardjfl and are adjustable in steps of about 
0.1 Hz and 0.02 degrees, respectively. The duration of each pulse is determined by 
a TTL output of the FPGA board with 10 ns resolution, which controls an rf switch, 
resulting in approximately square-shaped pulses. While the randomized sequences in the 
benchmarking experiment used at most 995 computational gates (1993 pulses), much 
longer sequences can be implemented with the present setup. 

The frequency of the qubit transition is acquired by deliberately detuning by about 
1 kHz, and performing a detuned Ramsey experiment (two 7r/2-pulses separated by a 
variable time). A fit to the resulting oscillation, with a delay time as long as 13.5 ms, 
yields the qubit transition frequency with statistical uncertainty less than 2 Hz. The 
measurement also yields an estimate of the unrefocused decay time T 2 * ~ 25 ms. 

We determine the pulse duration necessary for a tt/2 rotation by varying the time 
of the applied pulse and fitting to the measured Rabi oscillation. A typical tt/2 time is 
31.05 /xs, with an uncertainty of about 0.01 {is from the fit. 

We characterize the effect of systematic frequency control errors in the 
benchmarking sequence by performing 16 randomized computational gates sequences 
(4 CG and 4 PR) at a truncation length of 500 for various detunings. The CG and PR 
sequences used here are different than those used in the benchmarking protocol above. 
In Fig. [3](a) the average output state fidelity at each detuning is plotted. The output 
state fidelity is peaked at a detuning near zero (—10(7) Hz), confirming the accuracy of 
the measured qubit transition frequency. A gaussian fit yields a width of 150(13) Hz. 

d Incorporating an AD9951 DDS IC. 
e Using an Altera Cyclone II IC. 
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Figure 3. Evaluation of gate control errors using 16 randomized sequences truncated 
at 500 randomized computational gates, (a) Fidelity as a function of detuning from 
the measured |0) to |1) transition. A gaussian fit to the data gives a peak at —10(7) Hz 
and width 150(13) Hz. We estimate systematic detuning contributes < 1 x 10 -5 to the 
average error per gate, (b) The average fidelity as a function of the CG pulse duration, 
where each PR pulse is twice the CG pulse duration. Fitting to a gaussian yields a 
peak at —0.06(4) /is and a width of 1.1(1) /is. We thereby approximate that < 1 x 10 -5 
of the average error per gate results from systematic incorrect pulse duration. Error 
bars are the standard deviation of the mean. 

Given the uncertainty in the qubit transition frequency above and the deviation from 
zero of the fitted peak detuning, we estimate the contribution to the average error per 
gate resulting from systematic detuning as < 1 x 10 -5 . 

A similar evaluation is performed for the duration of the microwave pulses. Now 
we measure the average output state fidelity as a function of the pulse durations, with 
the result shown in Fig. [3](b). Again, the fidelity is peaked near the independently 
determined n/2 time of 31.05 /xs. Since the uncertainty in the nominal pulse duration 
was only 0.01 /xs, and a gaussian fit to the data in Fig. [3](b) gives a width of 1.1(1) /xs 
and a peak at —0.06(4) /xs, we expect systematic inaccuracies in pulse length contribute 
< 1 x 10 -5 to the average error per gate. 

The effect of pulse amplitude noise is estimated by measuring the power spectral 
density of continuous microwave amplifier output. Since fluctuations much faster than 
the pulse length will average to a consistent constant value, and systematic drifts with a 
time scale much longer than the pulse duration should be manifest in a non-exponential 
decay of the output state fidelity (absent in the data), we conservatively estimated the 
potential effect of amplitude noise in the range of about 1 kHz to 100 kHz. We estimate 
the total contribution to the average error per gate due to amplitude fluctuations in this 
range is < 1 x 10 -5 . 

Decoherence of the qubit is another error mechanism, which can arise as a result of 
inhomogeneity in the system or uncontrolled interactions such as spontaneous scattering 
of the optical lattice light. Dephasing due to static inhomogeneity can be counteracted 
by using a refocusing, or spin-echo, pulse. While the benchmarking protocol does not 
explicitly refocus the qubit, some refocusing is expected from the Pauli randomization 
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Figure 4. Contribution of decoherence to gate error, (a) Schematic depiction of 
determining the decoherence rate in the lattice using a standard detuned spin-echo 
technique, where for each time T the amplitude of oscillation is obtained by scanning 
St over about 1.5 ms. (b) Measured decay in the amplitude of the detuned spin-echo 
signal as a function of the time T, yielding an average decoherence time T 2 « 0.28(2) s 
(exponential fit), (c) Experiment schematic for measuring the contribution of the 
duration of the benchmarking sequence to the average error per gate. PR is a Pauli 
randomization (tt) pulse; CG is a computational gate (tt/2) pulse, (d) Average 
fidelity as a function of the time between pulses in the benchmarking sequence, th, 
at a truncation length of 500 randomized computational gates (1003 pulses). An 
exponential fit to the average fidelity versus total sequence time yields a decay constant 
of 0.31(4) s, consistent with T 2 . Uncertainties are the statistical standard error. 



pulses. 

We measure the decoherence time T 2 with a detuned spin-echo experiment, as 
illustrated in Fig. 0(a). The decay time is measured at two different magnitudes of 
the external magnetic field with similar results (average decay time T 2 ~ 0.28(2) s), 
indicating that the increased field sensitivity (from 4 kHz/mT to 37 kHz/mT) of the 
qubit states at the larger magnetic field does not limit the decoherence rate (Fig. BJb)). 
The data presented in Fig. [2] and [3] was taken at 0.32 mT. 

A more direct measure of decoherence in the presence of the partial refocusing from 
the benchmarking sequence is achieved by simply extending the total time over which 
the sequence of randomized computational gates is applied. As displayed in Fig. H](c), 
we accomplish this by adding an additional time th between subsequent pulses in the 
benchmarking sequence. Similar to the control error evaluation, a total of 16 randomized 
computational gate sequences (4 PR, 4 CG) at a truncation of 500 are applied to the 
atoms, with the average output fidelity for each value of th shown in Fig. S](d). A simple 
exponential fit to the average fidelity as a function of the total pulse sequence time 
gives a decay constant of 0.31(4) s, consistent with the spin-echo measured T 2 average 
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value of 0.28(2) s. Given that a single randomized computational gate (PR and CG) is 
typically 94 /xs, a T 2 decay during this time would result in an error per gate of about 
1.6(3) x 10 -4 . Thus, the primary contributor to the error in a single qubit gate appears 
to be the T 2 coherence time. 

A simple Monte Carlo simulation of the benchmarking sequences is used to model 
the effective refocusing of the qubit. Although the benchmarking sequences do not 
explicitly refocus the qubit, the simulation indicates that the qubit is almost completely 
refocused for a sufficiently long truncation of the benchmarking sequence. From the 
Monte Carlo simulation, we estimate the contribution of inhomogeneity across the atom 
cloud to the error per gate to be < 1 x 10 -5 , confirming that the dominate contribution 
to the error is the decoherence characterized by T 2 . 

5. Outlook 

Randomized benchmarking on atomic qubits in an optical lattice reinforces the viability 
of this system for applications in quantum information. A series of experiments 
demonstrates that we have precise control of the single qubit gate operations, without 
requiring composite pulse sequences or employing pulse shaping techniques. Although 
these methods have been shown to reduce the influence of control errors 0, M, CEE1 [EH fT8] 
and may be useful in future implementations, the gate operations are currently limited 
by the coherence time (T 2 ) of the system, which is probably constrained by the optical 
lattice. 

Spontaneous Raman scattering from the lattice is estimated from experimental 
parameters to occur at a rate of about 0.2 s _1 . Assuming that a single Raman scattering 
event depolarizes the qubit while Rayleigh scattering retains the coherence [TH [20] . 
spontaneous scattering then results in an error of about 1 x 10 -5 per randomized 
computational gate. The fidelity of the single qubit operations is probably limited 
by other aspects of the current optical lattice arrangement, such as time-dependent 
frequency shifts resulting from lattice beam intensity fluctuations. The scatter in 
Fig. [2](e) could be a result of these coherent errors. The operation of the quantum gate 
may be improved by better intensity stabilization of the lattice beams, or by operating 
the lattice in a regime that is insensitive to differential scalar light shifts [211 E2, E31 [24] . 

Advanced applications in quantum information will require individual qubit 
addressing. The microwave gate operations demonstrated here can be made to address 
single lattice sites by integrating high numerical aperture optics in the system to allow for 
tightly focused addressing beams [25]. In this configuration, the average error per gate 
resulting from both spontaneous scattering in the presence of the addressing beam and 
off-resonant microwave transitions in neighboring lattice sites can be < 1 x 10 -4 [25l l26]FI 
Thus, the combined addressing beam and optical lattice system should be able to utilize 
the high fidelity quantum gate operations presented here in a scalable architecture, 

f In this estimate we assume lattice beams of wavelength 1.064 /im, an addressing beam operating at 
789.5 nm with beam waist 1 /im, and gaussian-shaped microwave pulses. 
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facilitating more advanced quantum information protocols. 
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